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The two-layer square lattice quantum antiferromagnet with spins | shows a magnetic order- 
disorder transition at a critical ratio of the interplane to intraplane couplings. We investigate the 
dynamics of a single hole in a bilayer antiferromagnet described by a t — J Hamiltonian. To model the 
spin background we propose a ground-state wave function for the undoped system which covers both 
magnetic phases and includes transverse as well as longitudinal spin fluctuations. The photoemission 
spectrum is calculated using the spin-polaron picture for the whole range of the ratio of the magnetic 
couplings. This allows for the study of the hole dynamics of both sides of the magnetic order-disorder 
transition. For small interplane coupling we find a quasiparticle with properties known from the 
single-layer antiferromagnet, e.g., the dispersion minimum is at (±7r/2, ±7r/2). For large interplane 
coupling the hole dispersion is similar to that of a free fermion (with reduced bandwidth). The 
cross-over between these two scenarios occurs inside the antiferromagnetic phase which indicates 
that the hole dynamics is governed by the local environment of the hole. 



I. INTRODUCTION 



Since the discovery of high-temperature superconduc- 
tivity, doped antiferromagnets (AF) have been studied 
intensively. It is widely accepted that many properties of 
the superconducting cuprates are determined bv,the hole- 
doped CuC>2 planes. A number of experimentsu indicate 
that the cuprates are near a quantum-critical point of 
antiferromagnetic instability: The undoped materials are 
known to be antiferromagnetic Mott-Hubbard insulators, 
whereas hole doping destroys the antiferromagnetic long- 
range order (LRO) and leads to superconductivity. The 
investigation of the interplay between long-range mag- 
netic order and quantum disorder is therefore of great 
theoretical interest. 

A model system which shows a quantum transition 
between an ordered and a disordered magnetic phase is 
the S = | bilayer antiferromagnetffll Here each of the 
two planes is composed of a nearest-neighbor Heisenberg 
model with coupling constant Jii . The spins of corre- 
sponding sites of each layer are coupled antiferromag- 
netically with a coupling constant J±. In the limit of 
small J±/J\\ the model describes two weakly coupled 
AF planes. At T = this system possesses AF long- 
range order and gapless Goldstone excitations. In the 
opposite case of large Jj_/J||, pairs of spins interact- 
ing via Jj_ form spin singlets being weakly coupled by 
Jii. Then the spin excitations are gapped triplet modes, 
there is no magnetic LRO. At a critical ratio (J±/J\\) c 
a quantum transition between the two phases occurs 
which is believed to be of the 0(3) universality classou. 
The bilayer antiferromagnet has been studied by a num- 
ber of numerical and analytical techniques. Quantum 
Monte Carlo calculationscm and series expansionsQ yield 
an order-disorder transition point of (J±/J\\) c ~ 2.5. A 



similar result has also been obtained analytically using 
a diagrammatic approach to account Jjor the hard-core 
interaction between triplet excitationaJ. Bond-operator 
mean-field theory has-.recently been applied to the bi- 
layer Heisenberg AFli3 and gives a transition point of 
(</i/Jil)c ~ 2.3. Note that Schwinger boson mean- field 
theorya predicts a very large value of LLl/J\\) c ~ 4.5, 
and also self-consistent spin-wave theoryfcm, which yields 
(Jx/J\\)c ~ 4.3, fails to reproduce the numerical results. 
Chubukov and Morro have pointed out that this discrep- 
ancy is due to the neglect of longitudinal spin fluctuations 
in the conventional spin-wave approach. 

In this paper we discuss the bilayer antiferromagnet at 
zero temperature with hole doping as an additional de- 
gree of freedom. We consider the standard t — J model 
on a bilayer square lattice consisting of N sites per plane, 
so the total number of lattice sites is 2N. Each pair of 
corresponding sites in different planes is considered to 
form a rung, so we have N rungs. More precisely, the 
Hamiltonian reads 

H= -t x ^(4,<A 2 , CT + ff.c.) 

Yj (4m,*Cjm,*+ H.C.) 
(ij)ma 

i 

+ J\\ J]] (Sim ' Sjm — ) • (1) 

(ij)m 

Here and in the following i and j denote rungs of the 
bilayer lattice, (ij) denotes a summation over all pairs 
of nearest-neighbor rungs, m = 1,2 is the plane index, 
so im denotes a lattice site. Si m is the electronic spin 
operator and rii m the electron number operator at site 
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im. The electron operators c\ m a exclude double occu- 
pancies, cj m o . = c\ m (1 — n im ^ a ). We choose the z-axis 
along the direction of the rungs, i.e., perpendicular to the 
planes. 

First we want to comment the zero-temperature phase 
diagram of the doped bilayer t — J model which has 
to our knowledge not been systematically studied up to 
now. The following non-thermal control parameters can 
be considered: the ratio Ji/J||, the doping level S, and 
the relative hopping strength t/J. The doped single- 
layer antiferromagnet is known to exhibit a strong depen- 
dence of magnetic properties on the hole concentration 
5. With increasing hole concentration the staggered mag- 
netization decreases and vanishes at a critical hole con- 
centration S c p£ a_few percent where the system becomes 
paramagnetidliro. (This is .eapsistent with experiments 
on high-T c superconductorstirll3. ) On the other hand, 
in the undoped limit a large interplane coupling J± also 
destroys the antiferromagnetic LRO as discussed above. 
Thus it is likely that the system shows (i) AF LRO ac- 
companied by spontaneous symmetry breaking and the 
existence of massless Goldstone modes in the region of 
very small doping and small J±. Otherwise it is param- 
agnetic; here one has to distinguish between (ii) a spin- 
gapped phase which occurs at small 5, small hopping tu, 
and large interplane coupling, and (iii) a gapless phase at 
larger doping and small interplane coupling. The gapped 
phase (ii) is dominated by interplane singlet pairs (incom- 
pressible spin fluid) whereas the gapless phase (iii) shows 
only weak interplane coupling and behaves mainly like a 
single layer system at S > 6 C (compressible spin fluid). 
(Issues like pseudogap features and transport properties 
will not be adressed here.) In the present work we restrict 
ourselves to the zero-doping limit (one hole) . So the mag- 
netic (bulk) properties will not be affected by the doped 
hole, i.e., we probe the transition between the antiferro- 
magnetic (i) and the gapped (incompressible) paramag- 
netic phase (ii). 

The dynamical properties of holes in high-T c super- 
conductors have been subject to a lot of experimen- 
tal and theoretical work. Angle-resolved photoemiasion 
(ARPES) experiments at zercO or small dopingjiljU in- 
dicate a quasiparticle band with a small bandwidth pro- 
viding evidence for strong electronic correlations in the 
high-X^ compounds. A large number of numerical and 
analytical studies reveal that a single hole in an AF spin 
background has non-trivial properties: The spectral func- 
tion consists of a pronounced coherent peak at the bot- 
tom of the spectrum and a incoherent background at 
higher energies. The coherent peak can be associated 
with the motion of a dressed hole, ^._-a-iioie surrounded 
by spin defects ("spin polaron") yyMM The hole 
motion in bilayer antiferromagnets has been studied in a 
few papers and only for a small parameter range: Using 
self-cojnsistent Born approximation (SCBA) it has been 
shownca that small interplane coupling only weakly a-f-i 
fects the hole properties. Exact diagonalization studiesc3 



have indicated a larger effect of the interplane coupling 
at finite doping. However, the system sizes accessible to 
numerical methods (2 x 8 in Ref. ^) do not allow for the 
study of small hole densities and arbitrary momenta. 

In this paper we prefer an analytical approach to the 
one-hole problem which is based on the picture aLthp 
spin-bag quasiparticle (QP) or magnetic polaronE3~E3. 
The magnetic background for the hole motion will be 
modelled within a modified bond-operator representa- 
tion of the spins on each rung. In the ordered phase one 
type of triplet bosons condenses which will be included in 
the modified basis operators. So our excitation operators 
continuously interpolate between the triplet excitations 
of a singlet product state and the transverse and longitu- 
dinal excitations of a Neel-ordered state. The deviations 
of the spin background caused by tiw-hele-jnotion are 
described by a set of path operatorsHlI3ErEI The one- 
hole sp ectr al function will be evaluated using a cupulant 
versionE§E3 of Mori-Zwanzig projection techniqueEa. 

The paper is organized as follows: In Sec. II we pro- 
pose a ground-state wave function for the undoped bi- 
layer antiferromagnet. In Sec. Ill we develop the Hamil- 
tonian for the doped system in this new operator repre- 
sentation and discuss the various hole motion processes. 
Sec. IV focuses on the motion of a single hole in an 
otherwise half-filled bilayer antiferromagnet. The one- 
hole spectral function for this case corresponds directly 
to the result of an ARPES experiment on an undoped 
system. After a brief sketch of the projection technique 
we investigate the spectral properties in dependence of 
the magnetic background, especially when crossing the 
phase boundary between gapped paramagnet and anti- 
ferromagnet. The results show that the character of the 
one-hole spectrum is mainly determined by the local spin 
environment of the hole whereas the spectrum changes 
only weakly when crossing the order-disorder transition. 
A conclusion will close the paper. 



II. UNDOPED BILAYER ANTIFERROMAGNET 



To begin with, we consider the case of half-filling. Us- 
ing bond operators we derive a wave function which de- 
scribes the ground states in both the quantum disordered 
phase (at large J±) and the Neel phase (at small J±). 
The wave function will be constructed within a cumu- 
lant formalismt§Ej. It bases on an expansion around a 
product of rung states which are singlets in the disor- 
dered phase and sums of singlets and z-type triplets in 
the ordered phase. This wave function will be used in 
Sec. IV as "magnetic background" to calculate dynami- 
cal properties of a single hole. 
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A. Rung basis and Hamiltonian 

In the limit of vanishing intraplane coupling, Jn/Jj_ — > 
0, pairs of spins on each rung form a singlet ground state. 
The excitations are localized triplets with an energy gap 
J_i_. Switching on Jm leads to an interaction of triplets 

4l«) 

iio) 
4io) 



on neighboring rungs and to a dispersion of the triplet 
excitations. To describe this dipcrizcd phase we employ 
a bond operator representationE2l for the spins per rung. 
For each rung i containing two S = 1/2 spins Sji, Si2 
we introduce operators for creation of a singlet and three 
triplet states out of the vacuum |0): 

= ^=(1 U>-UT», 
= ^=(1 TT>-|U», 
= ^=(1 TT> + IU», 

= 4(iu)+ut)). (2) 



- il,T »24 

- — fc f c f 

V *1,T »2,t 

- — fc f c f 

- — fc f c f 



+ 4,i4a)io> 
+ 4,i4, T )io) 



Then the following representation holdsE3: 

^il,2 = 2^"' =S ^ ia 4 s i — i£uf3~i~t\pti~i) ■ (3) 

The algebra of the operators {si,t ix ,t iy ,t iz } has to be 
specified to reproduce the correct algebra for the spin op- 
erators. Here bond operators satisfying-the usual bosonic 
commutation relations will be chosentJ. In order to en- 
sure that the physical states are either singlets or triplets 
one has to impose the condition: 



s i s i ~l~ ^ ] t ia t ia 1. 



The singlet product ground state for J\\/J± 
written as 



»o = 



IKl 



(4) 

can be 
(5) 



We can now substitute the operator representation (|3j) 
into the Heisenberg Hamiltonian for the bilayer and ob- 
tain the following Hamiltonian. H . -r Hq + Hi in the bond 
(rung) operator representationEM3: 



Hq - Jj_ haha 



(ij)a 



( 44 44 tjahp )■ 



(6) 



{ij)af3 



The sum («j) runs over pairs of neighboring rungs of the 
lattice. The ground state of -ffo is given by the product 
state \(j)o) defined in eq. (JsJ) - 

Approaching the quantum critical point the triplet ex- 
citations will become gapless at momentum (%, 7r), see 



e.g. Refs. y,|6|. The magnetic quantum phase transition 
to an antiferromagnetically ordered state includes spon- 
taneous symmetry breaking and can be described via the 
condensation of triplet bosons in one particular direction 
(which has to be induced by an infinitesimal staggered 
field). Neglecting fluctuations a symmetry-broken state 
with a condensate of z-type triplet bosons can be written 
as 

\4> ) ~ 

exp(A5> iC ^4)|0 o ) = II(4 + Ae iQR *4)|0) (7) 



where the spontaneous staggered magnetization points 
in z-direction. The exponential term accounts for boson 
condensation, the condensation amplitude is given by A. 
A = means rotational symmetry in spin space (then 
\4>o) = 1 0o) ) 5 A ^ breaks this symmetry. A = ±1 trans- 
forms the singlet product state |</>o) into a Neel state. 
Q = (n, 7r) is the AF ordering wave vector. 

The basic idea for a proper description of the excita- 
tions of the product state |</>o) (@) is now to transform the 
basis states on each rung. We replace the basis operators 
Si for singlets and ti z for z-triplets by 
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VTT>? 
l 



(8) 



The new basis per rung consists of the operators 
{Si, Si, ti x , t iy } which still satisfy bosonic commutation 
relations. For A = this basis reproduces the usual bond- 
boson basis. For non-zero A the state sj|0) interpolates 
between the rung singlet and a (local) Neel state, the 



state zj\0) is orthogonal to sj|0). The product state 
(0) can now be written as 

i0o> = n4°>- 



(9) 
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For | A | = 1 this state is an antiferromagnetic Neel state. 
Its excitations are given by 4x*i> t\ y s i , and zjs^ Here, 

{t\ x ± it\ )s i create one spin flip in one of the planes and 
can be regarded as transverse spin fluctuations. The op- 
erator zjs i flips both spins on one rung which can be 
interpreted as longitudinal spin fluctuation. The basis 
parameter A has to be determined separately which will 
be discussed in the next subsection. Below it is shown 
that in the ordered phase A varies continuously between 



A and B denote the sublattices of the square lattice, i.e, 
giQRi _ f or j £ A or B, respectively. For A = 
the Hamiltonian (Jl^) is identical to derived with the 
usual bond operators. For |A| < 1 the ground state of 
Hq defined in (|l0|) is given by the modified product state 
\4>o) (||)- It can be seen that Hi in (|ic| ) contains creation, 
hopping, and conversion terms of the three types of exci- 
tations {zi,ti X ,ti y }. Contrary to the Hamiltonian Hi in 
(^) where the triplet excitations can only occur in pairs, 
here also single z-type excitations can be created and de- 
stroyed, e.g., by J±(zls i + s\z i ). The effect of these terms 
is directly related to the basis parameter A and will be 
used to determine its value. Furthermore it can be seen 
that all terms creating z-type excitations out of \cj)o) van- 
ish for Jl = 0, |A| = 1. This means that the ground state 
in the limit of vanishing intcrplanc coupling J_l does not 
contain single longitudinal spin fluctuations. They will 
become, however, important for describing the destruc- 
tion of the antiferromagnetism with increasing Jj_ . 

Note that z-type bosons formally occur together with 
transverse fluctuations even at J± = 0, |A| = 1. This 
happens if a transverse fluctuation created by Jii hits 



A = 1 for J± = and A = at the transition to the 
disordered phase. Note that the idea of a condensate for 
one type of bosons has also been employed by Chubukov 
and Morrcl who applied a modified version of non-linear 
spin-wave theory to the problem of the undoped bilayer 
antifcrromagnet . 

Using (H) we can reformulate the Hamiltonian for the 
(undoped) bilayer Heisenberg model in terms of the new 
basis operators {si,Zi 1 ti X , ti y }. Note that the splitting of 
the Hamiltonian is different from (H). 



(10) 



a second transverse fluctuation being already present at 
the same rung, in other words, if spins on correspond- 
ing sites in each of the planes are flipped independently 
by Jii . These processes have to be included in an ex- 
act treatment of the ground-state wave function; within 
the approximation of independent bosons (which will be 
described in Sec. II C) they are neglected. 



B. Approximation for the ground state 

To obtain an approximation for the g-tound-state wave 
function we adopt a cumulant approached which has been 
successfully applied to a variety of strongly correlated 
systems. One starts from a splitting of the Hamiltonian 
into H = Hq + Hi where eigenstates and eigenvalues of 
Hq are known. The ground state |^o) of" H can be con- 
structed from the unperturbed ground state \<po) of Hq 
by application of a so-called wave operator_Q which con- 
tains the effect of Hi. It has been showno that fl can 
be written in an exponential form, = e s , where S in- 
troduces fluctuations into \4>o}- The ground-state energy 



Hq — ji (t\ x t ix + t\ y t iy + - — ^ 2 (zj 



^ = J ^ E 7TA^ eiQRl + 

i 

E ( (4r*JiC + 4y*J»)( 



SiSj - z lZj )+H.c. 



(ieA,j<=B) 



-i V f tt it t- t- -tt tt t- t • + He 

n / / \ b ix b jx L iy b jy b ix b jy b iy b jx ~ 11 •' 



(i£A,j<=B) 



+ 



Jl 



2(1 


+ A2) (■ 




J\\ 


2(1 


+ A2) (■ 




J\\ 


2(1 


+ A 2 ) 2 




J\\ 


2(1 


+ A 2 ) 2 



E ( & 5 } ((1 - ^ 2 )¥^ - Witjx) + Ay^) ((1 - * 2 )¥« - 2\Zit jv ) + H.c. ) 
E ( *M (C 1 - x2 W jx + 2As i t J - x ) + t\ y ~z] ((1 - X 2 )z i tj y + 2As i t iy ) + H.c. ) 



(i£A,jEB) 



]T ((1 - \ 2 ) 2 {ztz)~ Si ~ Sj + zls]s iZj ) + 2A(1 - \ 2 ){s\z) ~ z\^ j )s^ J + H.c. ) 



(iGAJGB) 
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E is calculated^ according to E = ((j)o\HD,\(t>o) c where 
((f>o\...\4>o) c denotes a cumulant expectation value with 
respect to \4>o}- 

For the disordered phase of the system it is suitable 
to start from the singlet product state \(f>o) given in eq. 
(H) and to include fluctuations in form of pairs of triplet 
excitations. An appropriate ansatz for the wave function 
is 

l^o) = n I0o> 

= exp ^VijatiJjaSjSi + J2@ nSn \<M ■ ( n ) 
\ ija ™ / 

The first term in the exponential contains pairs of triplets 
(with arbitrary distance Rj — Rj ) . The operators S n in 
the second term represent higher-order fluctuations, e.g., 
four-triplet operators. They will be used in Sec. IV. For 
the determination of the coefficients /iy a and /3 n the cu- 
mulant method provides a set of non-linear equationscj, 



t> \(sls]t ja t ia y HQ\ 
b Q \SlHn\cbo) c - 



The dot • indicates that the quantity inside (...)' has to 
be treated as a single entity in the cumulant formation. 
The equations (|lj) are derived from the requirement that 
fi|0o) is an eigenstate of H. From symmetry it follows 
that the coefficients fiij a — fiij only depend on the dif- 
ference vector Hi — Rj (but not on a). So the state ji/'o) 
obeys rotational invariance in spin space. When evalu- 
ating the terms in ( |l2] ) the exponential series terminates 
after a few terms since the exponent contains only op- 
erators which create (but not destroy) fluctuations. The 
resulting non-linear coupled equations for fj,^ can be par- 
tially decoupled by a Fourier transformation to momen- 
tum space. The remaining set of equations has to be 
solved self-consistently. 

For the ordered phase the treatment is similar, how- 
ever, one has to take care of the broken rotational sym- 
metry. Starting point is the product wave function |</)o) 
(^) . In the operator £1 we now explicitely distinguish be- 
tween pairs of transverse (t\ x tj x , tl y ij ) and longitudinal 
(zjzj) fluctuations: 



(12) 



IV'o) = V \<h) 



exp 



f ft 4- ft ft 'Is .5. 



(13) 



The operators S n again contain higher-order fluctuations 
and will be discussed later. The parameters /xy , z/y , and 
n have to be determined in analogy to ( |l2| ) by the equa- 
tions 

0= ^ \{s\s)t jx t ix yHn\^ ) c , 

= ( fo\{s\s\z&y Hn\fo) c , 

o = (4> \slHn\4> ) c (w) 

which again reduce to (|l^) in the case of A = 0. For A = 
we expect the transverse and longitudinal fluctuations to 
become equivalent, fi^ — i/y. For |A| = 1 the amplitude 
of the longitudinal fluctuations in the ground state van- 
ishes, vij = 0, if one neglects higher-order processes in J» 
as noted above. Note that Q in ( [l3"| ) does not contain op- 
erators for the creation of single z-type excitations (zjsj) 
although such operators are present in the Hamiltonian 
(|j~0|). The reason is that such operators change the den- 
sity of the condensate of the z triplet bosons and refer 
to the same degree of freedom as the parameter A in the 
product state |0o)- The aim here is, however, to describe 
the condensation of the z-type bosons completely via the 
introduction of A and the transformed basis states. In 
this way Q, contains only multi-spin fluctuations. 

Now we discuss the parameter A which is expected to 
be zero throughout the disordered phase and to vary con- 
tinuously from to 1 with decreasing J± in the ordered 
phase. Within the cumulant method A can be determined 
consistently using the identity 



Q=(4>o\(4z l Y Hn\fo) c . (15) 

This can be understood as the condition that fl con- 
tains no additional term like £:j's i in the exponential, 
i.e., the boson condensation is fully accounted for by the 
transformed basis state \<fio). Within the cumulant for- 
malism this can be seen from the equality 

n\4> ) c = f>|e At Q^ ) c = tte At Q^ ) c , 

tq,=5> IQR *ti,. (16) 

i 

Here we have used the fact that a cumulant expectation 

value does not change when the operator e Atc * z is subject 
to cumulant ordering instead of being applied directly to 

the wave function (see Ref. ^). The operator f2e A *Q* 
can be interpreted as new wave operator (applied to |0o))- 
It now contains the additional parameter A which has to 

be determined from = (^oKsfSj' H fle xtc * z \(f>o) c which 
is equivalent to eq. (|l5|). 



C. Ground-state properties 

If we restrict the fluctuation operators in the expo- 
nential of £1 to excitation pairs as described above (i.e. 
neglect further operators S n ) and take fully into ac- 
count the constraint (Eh, we find a critical coupling of 
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(J±/J\\)c ~ 3.2 for the order-disorder transition. At this 
level of approximation the magnetization in the limit of 
decoupled planes (J± — > 0) is around 86 % of its clas- 
sical value being larger than predicted by Monte-Carlo 
and spin-wave calculations. The described approxima- 
tion can be systematically improved by including higher- 
order fluctuations. This behavior is also known from 
calculations withia. the coupled-cluster method for the 
Heisenberg modelo. It can be shown that in the limit of 
an infinite set of operators the cumulant method yields 
exact results. For instance, if we include two additional 
operators S n with three and four fluctuations on adjacent 
sites into the operator f2 (jl^) the transition point shifts 
to (J±/J\\)c ~ 2.9. However, the analytical and numer- 
ical effort to set-up and solve the non-linear equations 
equivalent to (|l^,|lj) increases drastically with including 
higher-order fluctuations. 



The solution of the non-linear equations ( jl2] , |l4| ) shows 
that the absolute values of all coefficients in f2 are small 
compared to 1. In the disordered phase the maximum is 
taken by Hio ~ 0.07 (the coefficient for nearest-neighbor 
triplet pairs) at the point of the phase transition. In 
the ordered phase the maximum is at J± = with 
/iio ~ 0.14. With increasing distance between the triplets 
the values decay fastly; higher-order fluctuations acquire 
a coefficient being one or more orders smaller in magni- 
tude. The fact that the coefficients are small is equivalent 
to the statement that the total triplet density even at the 
transition point is small (~ 0.1), see Ref. 0. 



In Fig. [j] we show results for the ground-state energy 
E and the staggered magnetization M for the approxi- 
mation using triplet pairs only as described above and for 
the spin- wave approach (see below). The data are com- 
pared with recent series expansions from Ref. [| It can be 
seen that all approximations well describe the behavior 
of the magnetization which first increases with increasing 
Jj_ and then drops to zero at the critical coupling. 



+ 

5; 
Si 




0.4 0.6 0.8 



1 

0.75 
0.5 
0.25 






b) : 






1 





0.2 



0.4 



0.6 



0.8 



FIG. 1. a) Re-scaled ground-state energy E and b) 
staggered magnetization M for different approximation lev- 
els of the ground-state wave functions ( |ll|) and (|l3|). Solid: 
Triplet pairs, no additional operators S n , constraint exact, 
(J±/J\\) c = 3.2. Dashed: Triplet pairs, constraint ignored 
(spin-wave), {J±JJ\\) C = 4. Dash-dot: Series expansion data 
taken from Ref. M for comparison, (J±/J«) c = 2.55. The ver- 
tical dashed line in a) indicates the position of the magnetic 
quantum phase transition. In b) the coupling dependence of 
the basis parameter A is also shown. 

Another approximation we want to discuss briefly is 
to neglect the constraint (|J) completelyotD and also the 
quartic triplet terms in the Hamitonians (^|) and (frC]). 
The parameter A is simply chosen so that the prefactors 
of the terms creating single z-bosons out of |</>o) vanish. 
One can formally set 5=1 (condensation of s-bosons); 
\<j>o) can be considered as "vacuum" . Neglecting the con- 
straint (||) implies that the three types of excitations 
{zi,ti x ,tiy\ are treated as independent bosons. After a 
Fourier transformation one arrives at 



jj(sw) _ 



X>£x*qx + 4t qy ) ((1 + A 2 )Ji + 8A 2 J,, + 2 7q (l - A 4 ) J,,) 



(1 + A 2 ) 2 

(1 + \ 2)2 E ~ z lK ((! ~ A ')^ + 16X2 J \\ + 2 7q(l - A 2 ) 2 J|| 



+ J ll E Tq( f q^-q- + *q»*-qi, + H C ' 



2 \2 



+ J\ 



(1-A 2 ) 
(1 + A 2 ) : 



E^q( 



4 5 -q + H -°- ) 



(17) 



G 



with 7 q = ^(cosq x + cosq y ) and 



A = 



/4J|| - J_l 
4Jn + Ji 



(18) 



in the ordered phase (4J|| > J±) and A = otherwise. 
The Hamiltonian ( |T"7|) can be easily diagonalized by a Bo- 
goliubov transformation. In this case the ground-state 
wave functions (11) and Jl3| ) (with pairs of triplet ex- 
citations) become exact. The described harmonic ap- 
proximation can be considered as linear spin-wave the- 
ory for the bilayer problem with inclusion of longitudi- 
nal spin fluctuations!! The critical coupling here becomes 
(J±/J\\) c = 4. In the limit of decoupled planes the results 
are equal to that of the linear spin- wave approach, e.g., 
the magnetization takes 60.5 % of its classical value, see 
Fig. 2. Further details on the ground-state calculations 
will be published elsewhere. 



III. DESCRIPTION OF HOLE MOTION 

If we remove one electron from rung i a single hole 
state on this rung is created. Let us denote by a\ m a the 
creation operator of a hole with spin a on rung i, where 
the index to = 1,2 represents the number of the plane. 
This operator is defined by the equation 



-J 

it n, a 



|o>, 



(19) 



where m = 3 — m. This is a fermionic operator which 
satisfies the usual anticommutation relations. The hole 
on the rung i interacts with the triplet excitations on its 
nearest-neighbor sites. The interaction Hamiltonian can 
be easily found by calculating all possible one-hole ma- 
trix elements of the initial Hamiltonian ([!]) (see Ref. 32). 
For the disordered phase (A = 0) the result is 



icr 

h \ - 
2 v 

(ij)mc 



+ H.c. 



a im,cr a jm,cr S j S i 



*n £ (* s,-, 

(ij) 



J\ 



+ H.c. 
H.c 



H S hj S j 



4 s, 



H.c. 



i Si 



t x t, 



H.c 



H.c 



(20) 



Following Rcfs. [32|,[36] we have introduced the notations 
t = (t x ,t y ,t z ) for the triplet vector and 



— ^ £ (~^) ma lm,a&aP a im,0, 

maj3 

— ^ £ a im,a^a0a jm l3 . 



(21) 



ma/3 



for generalized hopping operators where a is the vector 
of Pauli matrices. In addition we have to impose the con- 
straint that a hole state and one of the doubly occupied 
states cannot coexist on the same rung: 



sis. 



y ] " •m.n", 



= 1. 



(22) 



The several terms in the Hamiltonian ([20|) have been dis- 
cussed by EdeiE3. They contain hopping without chang- 
ing the spin background (direct hopping, second and 
third term) as well as spin-fluctuation assisted hopping 
(fourth and fifth term) and exchange processes where the 
hole remains on its rung (last two terms) . Together with 
the Hamiltonian for the doubly occupied rungs (|^) the 
Hamiltonian (^5|) contains the complete one-hole dynam- 
ics for the disordered phase. For the case of more than 
one hole additional hole interaction processes would have 
to be considered. However, in the present work we re- 
strict ourselves to the discussion of the one-hole problem. 

If we work in the magnetically ordered phase it is con- 
venient to use the generalized basis {s, z, t x , t y } from Sec. 
II to represent the states of the doubly occupied rungs. 
From the Hamiltonian ( |20| ) and the definition (||) of the 
basis operators one can construct a Hamiltonian for an 
ordered background state. It again contains processes 
proportional to t\\ where the hole hops to a neighboring 
rung and exchange terms propertional to Jy . To be short 
we only state some hopping matrix elements of the re- 
sulting Hamiltonian. For that we symbolically introduce 
two-rung states \XY) = XJyJ |0) where X and Y are the 
states on two adjacent rungs with i G A and j G B. 



(sa ]'\H h \a 01 s} 
(za ^\H h \a 01 z) 
(za ^\H h \a ^s) 
(t x aoi\H h \aois) 



h l - A 2 



2 1 



A 2 
A 2 



2 1 + A 2 ' 

*ll (1 + A) 2 
2 1 + A 2 
ill 1 + A 



2 v / TTA 2 



(23) 



The first two terms represent processes of hole motion 
without desturbing the spin background (so-called direct 
hopping) whereas the last two terms are examples for 
hole hopping with creation of a "spin defect". It can 
be seen that direct hopping is impossible in the limit of 
|A| = 1, i.e., in a Neel-ordered background. 
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IV. SINGLE HOLE DYNAMICS 

To investigate the hole motion we consider a one- 
particle Green's function describing the creation of a sin- 
gle hole with momentum k at zero temperature: 



gql,w) = (Mi 



i 



-Cko#o) 



(24) 



where z is the complex frequency variable, z = lo + irj, 
77 — ^ 0. The quantity L denotes the Liouville opera- 
tor defined by LA = [H, A]_ for arbitrary operators A; 
k = {k x ,k y ,k z ) is the hole momentum with k z = or 7r 
for the bonding or antibonding band, respectively. |^>o) is 
the full ground state of the undoped system as described 
in Sec. II. The effect of the hole creation operator Ci m ^ 
applied to |^o) can be related to the fermion operators 
a im a introduced in the last section. One obtains for ex- 
ample for i £ i sublattice 



5*i,tS!|0) - 
4 



1 + A 



\/2(l + A 2 ) 

^i,t4|o) = -^|4,tI ) 



4j°> 



(25) 



and analogous relations for other basis states and for 
i G B. 

The hole motion prqcesatjSL.will be described in the con- 
cept of path operators EjqEj which create strings of spin 
fluctuations attached to the hole. For the application of 
projection technique we define a set of path operators 
{Ai} which couple to a hole and create local spin defects 
with respect to the undoped ground state \ipo)- The first 
operator Aq is the unity operator, the second one A\ 
moves the hole by one lattice spacing creating one spin 
excitation and so on. We are interested in calculating dy- 
namical correlation functions for the operators {AjCka-}: 



G 



1 



la, J a' 



(z) = dj \(A lC ^ F (Ajc ka ,) |Vo> • (26) 



The Green's function ( |24| ) for the hole is then given by 
Goa,Qa{ z )- Using cunrulants the correlation functions G 
can be rewritten asEZI: 



G I<7 , Ja ,(z) = (0o|^(^/c ka )t 



1 



L 



n|0 o ) c . 

(27) 



The brackets (^o| ••• \4'o} c denote cumulant expectation 
values with \<j>o). The "unperturbed" ground state at 
half-filling, i.e., one of the product states (||) or @ for 
the disordered or ordered system, respectively, will be 
denoted in the following by |0o) f° r both phases. The 
operator f2 has been described above, it transforms \<f>o) 
being the ground state of Hq into the full ground state 
\4>o) oi H = Hq + Hi at half filling. In the following, f2 is 
approximated by an exponential ansatz as described in 



Sec. II. Since the expression (g7|) contains both f2 and 
fit, the expectation values can be calculated only in the 
lowest non-trivial order of the fluctuations in f2, i.e., the 
expansion will be restricted to linear fluctuation terms. 
In order to correctly describe the short-range spin fluctu- 
ations in the system, operators with up to 4 triplet excita- 
tions with a maximum distance of 4 lattice spacings have 
been employed. For the evaluation of their coefficients 
we have linearized the set of equations ( |L^ , |l~4| ). This is 
possible since the absolute values of the coefficients are 
small compared to 1 as discussed in Sec. II C. Therefore 
quadratic and higher terms in the fluctuations operators 
can be neglected. By comparing the coefficients of the 
short-range fluctuations obtained in this way with the 
ones from the full (non-linear) ground-state calculation 
of Sec. II C we have numerically verified that the error 
in the matrix elements ((3^) introduced by this lineariza- 
tion is smaller than a few percent. Since the linearized 
equations for the coefficients do not yield a phase transi- 
tion point we fix it at the value (J±/J\\) c = 2.55 known 
from series expansions^! and Monte Carlo calculations!]. 
(In the limit of an infinite set of operators {S n } within 
the cumulant method the solution of the non-linear equa- 
tions has to yield the "exact" transition point which we 
assume to be at (J±/J\\) c — 2.55.) In the ordered phase 
we set 



f (J±/J||)c- J±/J\\ 
(J±/J\\)c + J±/J\\ 



(28) 



which is suggested by the results obtained iu Sec. II. 

Using Mori-Zwanzig projection techniqueES one can de- 
rive a set of equations of motion for the dynamical cor- 
relation functions Gi a j a /(z). Neglecting the self-energy 
terms it reads: 



/J &Ko,I<7 (z) Gza,J<T' (z) = XKa,Ja> , 
la 

ttK*,J<T> {z) = zSkjSsc - U Ko,Lo" Xll», JcT > ■ (29) 



La" 



Xio-.ja' and LUja-^ja-i are the static correlation functions 
and frequency terms, respectively. They are given by the 
following cumulant expressions: 

Xla,Jv> = (M^( A lCka)' f {AjC k<7 >yn\(j) ) C , 

= faoin^jCfaO't (L(Ajc ka ,) )'r^ ) c . (30) 

These terms describe all dynamic processes within the 
subspace of the Liouville space spanned by the operators 
{AiCk a }. The use of cumulants ensures size-consistency, 
i.e., only spin fluctuations connected with the hole enter 
the final expressions for the one-hole correlation function. 

In the present calculations we have employed up to 
1600 projection variables with a maximum path length 
of 3. The neglect of the self-energy terms leads to a dis- 
crete set of poles for the Green's functions, so the present 
approach cannot account for lincwidths. In all figures we 
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have introduced an artificial linewidth of 0.2i|| to plot 
the spectra. For details of the calculational procedure 
see e.g. Rcf. p3. 



A. One-hole spectrum 

Now we turn to the discussion of the final results. 
First we consider the case of vanishing interplane hop- 
ping t± = 0, i.e., the hole motion is restricted to one 
plane. (In this case k 7 , can be dropped.) 




(7i/2,n/2) 



2-4 
Energy [f N ] 

FIG. 2. One-hole spectral function for t± = (no inter- 
plane hopping), a) J±/J\\ = 0, b) J±/J\\ = 1, and different 
momenta (k x ,k y ). The other parameters are J\\/t\\ = 0.4; the 
energies are measured in units of tu relative to the energy of 
a localized hole. 




FIG. 
J±/J\\ 



2-4 -2 

Energy [/„] 

3. Same as Fig. [| but for a) J±/J\\ = 2, 
= 2.4. 



The one-hole spectral function ImG(k,w) for t± = 0, 
J\\/tu = 0.4 and various values of J±/J\\ is shown in Figs. 
H - Q For zero and small Jj_, i.e., in the limit of decou- 
pled planes, we observe a pronounced quasiparticle (QP) 
peak and a weak background at higher energies. The QP 
peak follows a dispersion with minima at (±7r/2, ±7r/2). 
These features are well-known from the one-hole problem 
in the single-layer antiferromagnet. With increasing J± 



the dispersion is " washed-out" because the antiferromag- 
netic correlations in the background state are weakened. 
At J±/J\\ ~ 1.5 a cross-over to a dispersion with min- 
ima at (tt, tt) and (0, 0) occurs. Note that this cross-over 
point still lies inside the antiferromagnet phase. Near 
the phase transition the spectral weight of the lowest pole 
at (0, 0) becomes small. Most of the spectral weight can 
be found in a band with minimum at (tt, tt) and maxi- 
mum at (0,0). The weak band visible at the bottom of 
the spectrum around momentum (0, 0) can be considered 
as "shadow band" originating from the antiferromagnetic 
background, i.e., it is obtained by shifting the band with 
large weight by Q = (tt, tt). 




2-4 
Energy [{„] 

FIG. 4. Same as Fig. |, but for a) J±/J\\ = 2.7, b) 
Jx/J\\ = 10. 

Further increasing J_l drives the system into the 
gapped paramagnetic (PM) phase. At the phase tran- 
sition it can be seen that shadow bands being present 
in the AF phase disappear in the PM phase. At large 
J_l the hole behaves like a free fermion with a dispersion 
proportional to t\\ (cos k x + cos k y ). 

Now we consider the effect of interplane hopping. We 
fix the ratio of the parameters to 



Jj 



_!L 

J\\ 



u 
1 



(31) 



which follows from the derivation of the t — J model 
from a Hubbard model for the bilayer system with on- 
site repulsion U. Spectra for different parameter sets 
with U/U — 10 are shown in Figs. ||- 0- An interplane 
coupling of t±/tu = 1 has already changed the character 
of the bands compared to t\ — (cf . Fig. [||a) . Further 
increasing t±/t^ to 1.5 (Fig. ^) forms a pronounced band 
with a free- fermion form at the bottom of the spectrum in 
the bonding channel (k z = 0). In contrast, the antibond- 
ing spectrum (k z = 7r) becomes incoherent since spectral 
weight is transferred to higher energies. Since we are still 
in the antiferromagnetic phase the lowest peaks in the an- 
tibonding channel appear at the same energies as in the 
bonding channel shifted by momentum Q = (7r, 7r,7r), 
i.e., these can again be considered as shadow bands. At 
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t±/tu > y/2.55 the AF order is destroyed. As shown 
in Fig. for t±/tn = 2 the shadow bands are disap- 
peared, and in the antibonding channel a pronounced 
free-fcrmion band at high energies is formed. Higher val- 
ues of t j_ will completely suppress the incoherent weight 
left in the antibonding channel. 



(l(,0) 



(7C/2.7I/2) 



(0,n) 



(71/2,11) 



(71,71) 



(7t/2,n/2) 
(0,0) 

2 -4 
Energy [%] 

FIG. 5. One-hole spectral function for t±/tu = 1, 
J±/J\\ = 1) J|| A|| — 0.4, and different momenta. The left 
and right panel show the bonding (k z = 0) and antibonding 
(k z = 7r) bands, respectively. 
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-4 -2 2 4 -4 -2 2 
Energy [%] 

FIG. 6. Same as Fig. |, but for t±/t\\ = 1.5, 
J±/J|| = 2.25, and J||A|| = 0.4. 



Varying the ratio J\\/t\\ (which is 0.4 for all figures 
presented here) does not change the picture qualitatively. 
Larger values of Jh suppress the incoherent background in 
all spectra. For very small values of Jn the spectral func- 
tions become incoherent since the hole is dressed with an 
increasing number of spin fluctuations, i.e., the radius of 
the spin polaron increases. Furthermore, the bandwidth 
of the quasiparticle dispersion at small J± is mainly con- 
trolled by J|| (instead of tu). This means that the band- 
width in units of t\\ decreases with decreasing J\\/t\\, see 
also next subsection. The location of the cross-over be- 
tween the two dispersion forms depends weakly on J\\/t\\ 
which will be discussed below. 
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FIG. 7. Same as Fig. |, but for t±/t\\ = 2, Jx/J\\ = 4, 
and Ji|A|| = 0.4. 



B. Quasiparticle bands 

Next we are going to examine the properties of the 
low-lying bands. Again we first discuss the case of van- 
ishing interplane hopping, t±_ = 0. For J± <C Jn one 
finds a narrow QP band shown in Fig. ^ (top panel). It 
has minima at (±7r/2, ±7r/2) and a bandwidth of approx- 
imately 1.5Jy. It corresponds to the coherent motion of 
a dressed hole, i.e., a spin-bag quasiparticle. (Note that 
a larger set of path operators would give a larger band- 
width of around 2 Jn which is knowjiirom the calculations 
for a single-layer antiferromagneta. However, paths of 
length 4 and more are hard to access for the bilayer sys- 
tem due to the increasing numerical effort.) 

The main contribution to the hole motion in the AF 
case can be understood as follows: the hopping hole lo- 
cally destroys the antiferromagnetic spin order leaving 
behind a string of spin defects. Quantum spin fluctua- 
tions can repair pairs of frustrated spins which leads to a 
coherent motion in one of the antiferromagnetic sublat- 
tices (spin- fluctuation-assisted hopping). For tii/Ju > 1 
the bandwidth of this coherent hole motion is of order Jy 
because the spin-flip part of the intraplane Heisenberg 
interaction Jii is necessary to remove the spin defects 
caused by hopping. 

With increasing J± the bandwidth of the dispersion 
first decreases; at a value of J±/Jh ~ 1.5 the minima 
move to (7T, 7r) and (0,0). Slightly larger values of J± 
lead to an increase of the bandwidth, and spectral weight 
is transferred to a band which has nearly tight-binding 
form. So the lowest pole at J±/J\\ — 2.4 has very small 
weight around momentum (0, 0). This shadow band is 
shown as dashed line in Fig. H. In the paramagnetic 
phase for larger J± the shadow band disappears. The 
QP peak has a dispersion being nearly the one of a free 
fermion but with a reduced bandwidth. This behavior 
follows from the existence of the direct hopping term with 
prefactor tu/2 in Hh (20). It leads to a disperison pro- 
portional to the saturation bandwidth equals half of 
the bandwidth of the uncorrelated system (8i||). 
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(0,0) (71/2,11/2) (3l.1t) (n/2,7t) (0.71) (7t/2,7l/2) (TC.O) 
In-plane momentum 

FIG. 8. Quasiparticle dispersion for t± = (no interplane 
hopping), J\\/t\\ = 0.4, and different values of Jx/Jy. The 
energy zero level has been set at the center of mass of the 
band. For J±/J\\ = 2.4 the heavy line shows the dispersion of 
the peak carrying the main spectral weight whereas the thin 
line corresponds to the lowest pole in the spectrum (shadow 
band, see text). 




(0,0) (71/2,71/2) (71,71) (71/2,71) (0.71) (7t/2,7l/2) (71,0) 

In-plane momentum 



t^i t,| jj 

FIG. 9. Quasiparticle dispersion for -jj- = -j- — -j with 

U/t\\ — 10 and different values of . The in-plane momen- 
tum (k x ,k y ) varies along the horizontal axis. Solid/dashed 
lines correspond to bonding/antibonding modes (k z — or 

77). 

The results for non-zero interplane hopping t±_ are 



shown in Fig. ^. Again, for small t± (small J±) the 
dispersion character is " antiferromagnetic" whereas for 
large t± (large J±) the bands have a simple tight-binding 
form. For Jj_ Ju and ij_ their dispersion is given 
by 

£k = — 1± cos k z + £|| (cos k x + cos k y ) (32) 

The bandwidth for the in-plane motion is reduced by a 
factor of 2 compared to a free fermion since the relevant 
hopping matrix element is (sao^\Hh\ao^s) = tn/2. 




FIG. 10. Bandwidth of the quasiparticle dispersion 
(shown in Fig. ^) vs. J±/J\\ for t± = (no interplane 
hopping) and J\\/t» = 0.4. The vertical dashed line again 
indicates the position of the magnetic phase transition. 

In Fig. [l0| we finally show the bandwidth of the QP 
dispersion for the case of vanishing interplane hopping 
and = 0.4, i.e., for the dispersions shown in Fig. 

^ It can clearly be seen that the cross-over from an " an- 
tiferromagnetic" to a simple tight-binding dispersion oc- 
curs at J±/ 'J|| ~ 1.5 and is connected with a pronounced 
minimum in the QP bandwidth. For J±/J\\ > 1.5 the 
bandwidth increases and saturates for large J± at 4in 
as discussed above. Results for other values of J\\/t\\ arc 
qualitatively similar. The minimum in the dispersion (i.e. 
the cross-over point) shifts with J\\/t\\: for = 0.1 it 

is located at J±/Jn ~ 1.1, for J\\/t\\ — 10 it is found at 
J±/J\\ ~ 2.0. This in turn means that the cross-over in 
the shape of the single-hole dispersion can be driven by 
varying J\\/t\\ at fixed J±/J\\. 



C. Relation to the short-range spin correlations 

In order to understand the behavior of the QP dis- 
persion and its bandwidth we illustrate in more detail 
the connection between the in-plane hole motion and 
the short-range in-plane spin correlations. The basic in- 
gredience for the observed behavior of the QP disper- 
sion is the competition between direct hopping which 
dominates in the disordered phase and spin-fluctuation- 
assisted hopping known from the single-layer AF. The 
relevant matrix element for direct hopping between states 
without spin deviations is t\ := (ipolc^HtCialipo) where 
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i,j are nearest-neighbor sites and \ipo) denotes theuup. 
doped background state. This term can be expressecoS 
by the static in-plane nearest-neighbor spin correlation 
Sr = (^o|S • S R |^o) with R = (1,0) as 

ti = 2t||(Sio + i). (33) 

Spin-fluctuation-assisted hopping is more complicated 
since it involves two hopping steps and one spin- 
fluctuation process. The matrix element for the 
spin-fluctuation process (being the most important 
for not too small values of t\i) can be written as 
h ■= (^o\c] a HjA2C. la \ipo) where A 2 moves the hole by 
two hopping steps (creating a path of spin defects) and 
i,j are now next-nearest neighbors. Transforming the 
expectatiatt value into spin correlation functions one 
obtamsE2rc3 

+ + ^ (34) 

where the average over the possible paths of length 2 has 
already been performed. 





1 1 1 1 '. 1 1 1 1 

1 h = 1 | 


yT. ^10 


V n / = 0.4^ ■ 



01234012345 



FIG. 11. Left: J±/J\\ -dependence of the equal-time 
short-range spin correlation functions Sr in the spin back- 
ground state calculated with the linearized exponential ansatz 
used for the description of the hole motion. Right: A vs. 
J±/J\\ from eq. (g) for J\\/t\\ = 0.4 and 1.0. 

The in-plane dispersion shape can be described by 
the energy difference A = E(tt,-k) — E(n/2,7r/2). Val- 
ues A > correspond to an " antiferromagnetic" disper- 
sion whereas A < occurs for a nearest-neighbor tight- 
binding dispersion. Neglecting longer paths and more 
complicated contributions to the hole motion the QP dis- 
persion is given by the sum of a nearest-neighbor and 
a next-nearest neighbor dispersion originating from the 
two processes described above. Then A can be roughly 
estimated from the matrix elements t\ and t-2.'- 

A ~ -Ah + t 2 . (35) 

The prefactor of t\ arises from the number of nearest 
neighbors; for the ^-prefactor the influence of two hop- 
ping steps has to be kept in mind, a fit to numerical re- 
sults for the single-layer problem yields a value of order 1. 
From this we see that the cross-over phenomenon can be 



understood in terms of the interplay between short-range 
spin correlations and the ratio of J\\/t\\- 

Fig. [ll] shows the values of Sr. for nearest and next- 
nearest neighbor sites obtained from the present calcula- 
tion. With increasing J±/J\\ the absolute values of the 
in-plane spin correlation functions decrease from their 
single-layer values and are weakened within the antifer- 
romagnetic phase. At the transition point the magnitude 
has dropped by a factor around 2 — 3 compared to J± = 
which again coincides with the fact that the density of 
spin (triplet) excitations is small in the disordered phase 
even at the transition pointfl. We have also plotted the 
quantity A estimated from eq. ( |35| ) for two different val- 
ues of J\\/t\\- One notices a semi-quantitative agreement 
of the zero in A and the dispersion minimum described 
in Sec. IV B. 



V. CONCLUSIONS 

In this paper we have presented for the first time a 
systematic analytical study of the one-hole dynamics on 
both sides of a magnetic ordering transition in a low- 
dimensional antiferromagnet. The system under consid- 
eration was a bilayer antiferromagnet described by a t- J 
Hamiltonian in the limit of zero doping. The magnetic 
background state has been modelled with modified bond 
operators. In the disordered phase these operators de- 
scribe the singlet ground state and the triplet excita- 
tions. In the AF ordered phase they account for the 
condensation of one type of triplet bosons and describe 
transverse as well as longitudinal fluctuations. Using the 
spin-polaron concept which describes the one-hole states 
in terms of local spin deviations we have calculated the 
one-hole spectral function for the whole range of mag- 
netic couplings. 

For the disordered background (large J±/J\\, gapped 
spin excitations) spin fluctuations around the hole are 
suppressed. The hole motion is dominated by direct 
hopping processes, i.e., hopping without disturbing the 
spin background. In the ordered phase for very small 
interplane coupling we recover the results known from 
the single-layer hole motion: The spectrum consists of 
a coherent QP peak at the bottom and an incoherent 
background. The QP can be associated with a mobile 
hole dressed by spin fluctuations. The bandwidth of its 
coherent motion is controlled by Jii . 

The cross-over between these two scenarios occurs 
inside the ordered phase where the antiferromagnetic 
short-range correlations become weakened. The cross- 
over is located between 1 < Jj_/J\\ < 2 depending on 
J||/t|| (for J||/i|| = 0.4 it is found at J±/J\\ ~ 1.5). 
Note that the cross-over can also be driven by varia- 
tion of the hopping strength i|| at a fixed value of J±/J\\ 
well in the antiferromagnetic phase (e.g. 1.3). This 
behavior follows from the competition between direct 
nearest-neighbor hopping and spin-fluctuation-assisted 
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next- nearest-neighbor hopping. In contrast, when cross- 
ing the magnetic phase boundary at J±/J\\ ~ 2.5 there 
are no drastic changes in the spectrum (and therefore in 
the ARPES response of such a system) . The only differ- 
ences between the spectra in both phases near the phase 
transition are weak shadow bands in some regions of the 
Brillouin zone in the AF phase. 

Note that our approximation which describes the hole 
motion processes in terms of short-range spin fluctua- 
tions is questionable in a region in the close vicinity of the 
phase transition due to the existence of long-ranged crit-. 
ical fluctuations. However, it has been argued recentlyE3 
that the influence of the critical modes suppresses the 
quasiparticle weight only in the limit of vanishing hop- 
ping t/J — ► 0. In contrast, for finite hopping the number 
of spin fluctuations near the hole remains finite even at 
the transition point. So we expect the picture presented 
in this paper to be valid at least in the regions away from 
the transition, i.e., \ J±/J\i — 2.55 1 > 0.1, which is sup- 
ported by the fact that our spectral functions on both 
sides of the phase transition (e.g. at J±/J\\ = 2.4 and 
2.7) do not show major differences. 

So the main result of the present paper can be sum- 
marized as follows: A magnetic phase transition has only 
weak influence on the ARPES spectrum of a doped an- 
tiferromagnet. The properties of the spectrum are, how- 
ever, dominated by the short-range environment of the 
hole. The statement is expected to hold also for the 
doping-induced phase transition in the single-layer AF. 
Here the antiferromagnetic correlations are strong even 
in the paramagnetic phase, i.e., from the above consid- 
erations one expects to find an "antiferromagnetic" hole 
dispersion on both sides of the transition. . This is ex- 
actly what was observed in recent workciHi3: analytical 
investigations of one hole in the AF phase as well as nu- 
merical studies of finite systems (which have a singlet 
ground state without long-range order) both show a co- 
herent hole motion with a dispersion of width 2J and 
minima at (±7r/2, ±7r/2). 

It should be pointed out that the above discussion ap- 
plies to intermediate and high energy scales (order t, J) 
only. Of course there exist low-energy properties of the 
spectrum which are expected to be influenced by quan- 
tum criticality, e.g., the lincwidth in a finite-temperature 
photoemission experiment should show scaling behav- 
ior in thc_nuantum-critical region associated with the 
transitiorO. These features as well as the hole dynamics 
in the bilayer system at low but finite doping are be- 
yond the scope of the present study and will be subject 
of future research. 
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